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' (— 1 ' In view of the latest T2K and MINOS neutrino oscillation data which hint at a rela- 

tively large 13 , we perform a systematic study of the Majorana neutrino mass matrix M v 
Qh| with two independent texture zeros. We show that three neutrino masses (m 1 ,m 2 ,m 3 ) 

^ e— I and three CP-violating phases (5, p, a) can fully be determined from two neutrino mass- 

squared differences (5m 2 , Am 2 ) and three flavor mixing angles (#12) ^23) $13)- We find 
that seven patterns of M v (i.e., A 12 , B x 2 3 4 and C) are compatible with current ex- 

■^j- ■ perimental data at the 3<r level, but the parameter space of each pattern is more strictly 

constrained than before. We demonstrate that the texture zeros of M v are stable against 

■^j - ■ the one-loop quantum corrections, and there exists a permutation symmetry between 

Patterns A 1 and A 2 , B^^ and B 2 or B 3 and B 4 . Phenomenological implications of M v 
on the neutrinoless double-beta decay and leptonic CP violation are discussed, and a 
1 . realization of those texture zeros by means of the Z n flavor symmetries is illustrated. 
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1 Introduction 



Compelling evidence in favor of neutrino oscillations has been accumulated from a number 
of solar, atmospheric, reactor and accelerator neutrino experiments since 1998 [1]. We are 
now convinced that the three known neutrinos have nonzero and nondegenerate masses, and 
their flavor states can convert from one kind to another. In particular, two neutrino mass- 
squared differences (5m 2 , Am 2 ) and two flavor mixing angles (0 12 , 9 23 ) have been determined 
to a reasonably good degree of accuracy from currently available experimental data. In spite 
of this progress made in neutrino physics, our quantitative knowledge about the properties of 
massive neutrinos remain quite incomplete — for instance, the absolute mass scale of three 
neutrinos, the sign of Am 2 and the smallest flavor mixing angle 9 13 are still unknown. Although 
the nature of massive neutrinos is also an open question, we assume them to be the Majorana 
particles. In this case the 3x3 Majorana neutrino mass matrix M v is symmetric, and the 
corresponding 3x3 flavor mixing matrix V contains three CP- violating phases (5, p, a) which 
are entirely unrestricted by current experimental data. 

The T2K [2] and MINOS [3] accelerator neutrino experiments have recently announced 
their preliminary data on the appearance v — > u e oscillation, which hints at a relatively large 
value of 6> 13 . If this mixing angle is confirmed to be not very small by new data from these two 
experiments and from the upcoming reactor antineutrino experiments (e.g., the Double Chooz 
experiment in France [4] , the Daya Bay experiment in China [5] and the RENO experiment in 
Korea [6]), it will be a great news to the long-baseline neutrino oscillation experiments which 
aim at a determination of the sign of Am 2 and the magnitude of the Dirac CP-violating phase 
5. The observed values of relevant neutrino flavor parameters may then allow us to reconstruct 
the Majorana neutrino mass matrix M v in the flavor basis where the charged- lepton mass 
matrix M l is diagonal: M v = VM V V T , where M v = Diag{m 1 , m 2 , m 3 } with m i (for i — 1, 2, 3) 
being the neutrino masses, and V is just the 3x3 neutrino mixing matrix. 

In the lack of a convincing flavor theory, at least five approaches have so far been tried to 
study the flavor problems of massive neutrinos [7] : (a) radiative mechanisms, (b) texture zeros, 
(c) flavor symmetries, (d) seesaw mechanisms, and (e) extra dimensions. Some of them are 
certainly correlative. For example, the neutrino mass matrix M u may have a few texture zeros 
as a natural consequence of an underlying flavor symmetry in a given model with or without 
the seesaw mechanism. Such texture zeros are phenomenologically useful in the sense that 
they guarantee the calculability of M v from which both the neutrino mass spectrum and the 
flavor mixing pattern can more or less be predicted. Note that texture zeros of a fermion mass 
matrix dynamically mean that the corresponding matrix elements are sufficiently suppressed 
in comparison with their neighboring counterparts [8], and they can help us to establish some 
simple and testable relations between flavor mixing angles and fermion mass ratios — if such 
relations turn out to be favored by the experimental data, they might have a fundamental 
reason and should originate from the underlying flavor theory. Hence a phenomenological 
study of possible texture zeros of M v does make a lot of sense. 

As M u is symmetric, it has six independent complex entries. If n of them are taken to be 
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vanishing (i.e., M v has n independent texture zeros), then we shall arrive at 

n n\(6-n)\ 



different textures. It is easy to show that a texture of M v with more than two independent 
zeros (i.e., n > 3) is definitely incompatible with current experimental data on neutrino masses 
and flavor mixing angles [9]. Hence a number of authors have paid particular interest to the 
two-zero textures of M u [10] — [14] 1 . There are totally fifteen two-zero textures of M v , which 
can be classified into six categories: 
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in which each "x" denotes a nonzero matrix element. Previous analyses of these fifteen patterns 
(see, e.g., Ref. [13]) led us to the following conclusions: (1) seven patterns (i.e., A 12 , B x 2 3 4 



: It is worth mentioning that the one-zero textures of M v have much less predictability than the two-zero 
textures of M v , and their phenomenological implications have been discussed in Ref. [15]. 
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and C) were phenomenologically favored; (2) two patterns (i.e., D x 2 ) were only marginally 
allowed; and (3) six patterns (i.e., E x 2 3 and F t 2 3 ) were ruled out. Taking account of the latest 
T2K and MINOS neutrino oscillation data, one may wonder whether the above conclusions 
remain true or not. A fast check tells us that Patterns D x 2 can be ruled out by today's 
experimental data at the 3a level, simply because 9 12 > 38° is no more favored. In addition, 
Patterns E 12 3 and F 1 2 3 remain phenomenologically disfavored. The question turns out to be 
whether Patterns A 12 , B 1)2) 3 )4 and C can survive current experimental tests and coincide with 
a relatively large value of 9 13 . The main purpose of this paper is just to answer this question. 

We aim to perform a very systematic study of the aforementioned seven patterns of M v 
with two independent texture zeros based on a global fit of current neutrino oscillation data 
done by Fogli et al [16] 2 . Because a two- zero texture of M v is very simple and can reveal the 
salient phenomenological features of a given Majorana neutrino mass matrix, we intend to take 
such a systematic analysis as a good example to show how to test possible textures of M v by 
confronting them with more and more accurate experimental data. Hence this work is different 
from any of the previous ones. In particular, almost all the physical consequences of Patterns 
A 12 , B 1 234 and C are explored in an analytical way in the present work; the stability of 
texture zeros and the permutation symmetry between two similar patterns are discussed; the 
numerical analysis is most updated and more complete; and a realization of texture zeros by 
mean of certain flavor symmetries is illustrated. We find that Patterns A 12 , B 12)3 4 and C 
of M v can all survive current experimental tests at the 3a level, but we believe that some of 
them are likely to be ruled out in the near future. 

The remaining parts of this paper are organized as follows. In section 2 we show that the full 
neutrino mass spectrum and two Majorana CP- violating phases (p, a) can all be determined 
in terms of three flavor mixing angles (0 12 , 23 , 13 ) and the Dirac CP-violating phase (5) for a 
given two- zero texture of M v . Section 3 is devoted to a complete analytical analysis of Patterns 
A 1 2 , B x 2 3 4 and C of M u , and to some discussions about the stability of texture zeros against 
the one-loop quantum corrections and the permutation symmetry between two similar patterns 
under consideration. With the help of current neutrino oscillation data we perform a detailed 
numerical analysis of those typical patterns of M v in section 4. Section 5 is devoted to some 
discussions about how to realize texture zeros of M v by means of certain flavor symmetries, 
and section 6 is a summary of this work together with some concluding remarks. 

2 Basic Formulas 
2.1 Important relations 

In the flavor basis where the charged-lepton mass matrix M t is diagonal, the Majorana neutrino 
mass matrix M v can be reconstructed in terms of three neutrino masses (m l ,m 2 ,m 3 ) and the 

2 Note that Schwetz et al have done another global fit of current data [17]. Although their best-fit results of 
three flavor mixing angles are slightly different from those obtained by Fogli et al [16], such differences become 
insignificant at the 3a level. Therefore, we shall mainly use the 3a results of [16] in our numerical calculations. 
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flavor mixing matrix V. Namely, 



M=V 



m 2 




V 1 



(8) 



It is convenient to express V as V = UP, where U denotes a 3 x 3 unitary matrix consisting 
of three flavor mixing angles (0 12 , 9 23 , 9 13 ) an d one Dirac CP-violating phase (5), and P = 
Diag{e v , e la , 1} is a diagonal phase matrix containing two Majorana CP- violating phases (p, a). 
More explicitly, we adopt the parametrization 



U = 



( c 12 c 13 

— C 12 S 23 S 13 — S 12 C 23 e 

\ — CinCnoSio + SloSnqe 



s 12 c 13 
— • S 12 S 23 S 13 + C 12 C 23 e 



" S 12 C 23 S 13 



C 12 S 23 e 



Coo C 



13 \ 

3 



S 23 C 13 



(9) 



23^13 



where s„- 



sin 9, :A and c,„- = cos 9 {j (for ij = 12, 23, 13) are defined. The neutrino mass matrix 



M v can equivalently be written as 



M„ 




(10) 



where X 1 = m 1 e 2lp , A 2 



m e 



2ia 



and A 3 = m 3 . If two independent elements of M v are vanishing 



(i.e., (M u ) ab = (M u ) a/3 = with ab 7^ af3) as shown in Eqs. (2) — (7), one can obtain [11] 

Aj ^3^63^2^/32 ~ ^2^)2^3^/33 

^ UaUnUrtUp - U al U bl U a2 U 02 ' 

A 2 ^al^fcl^n3^/33 ~ ^a3^fc3^nl ^/31 

from which two neutrino mass ratios are given by 



(11) 



m, 



mo 



/33 



c = — 

m 3 



U al U bl U a3 Up 3 - u a3 u b3 u al u pi 



U^U^UaiUfn — U al U bl U a2 Up 2 
and two Majorana CP-violating phases turn out to be 

^o3^b3^a2^g2 ~ ^a2^b2^a3^/33 

U^U^U^U^ — U al U bl U a2 U p 2 

'u al u bl u a3 u P3 -u a3 u b3 u al u pi 



(12) 



p = 



(7 = 



arg 
arg 



L u a2 u b2 u al u fil - u al u bl u a2 u t 



02 



(13) 



2.2 Parameter counting 

Since we have assumed massive neutrinos to be the Majorana particles, there are nine physical 
parameters: three neutrino masses (m 1 ,m 2 ,m 3 ), three flavor mixing angles (#125 #235 $13)5 an d 
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three CP- violating phases (5, p, a). By imposing two independent texture zeros on M v , we 
obtain four constraint relations as given in Eqs. (12) and (13). With the help of current 
experimental data on three flavor mixing angles and two independent neutrino mass-squared 
differences, defined as [16] 

5m 2 = m\ — m\, Am 2 = m 2 — - (m\ + m 2 ) , (14) 

one may determine or constrain both the neutrino mass spectrum and three CP-violating 
phases through Eqs. (12) and (13). 

We take (9 12 , 9 23 , 9 13 ) and (5m 2 , Am 2 ) as observables to see why the neutrino mass spectrum 
and three CP-violating phases can in principle be determined or constrained. First of all, note 
that £ = m 1 /m 3 and ( = m 2 /m 3 are functions of the Dirac CP- violating phase 5 as shown in 
Eq. (12). Then it is possible to determine or constrain 5 from the relation 

R - 5m2 2 (c 2 -a (15) 

u ~\Am 2 \ |2-(C 2 + £ 2 )|' 
Once 5 is fixed, we can obtain (£, () and (p, a) from Eq. (12) and Eq. (13), respectively. In 
addition we have 

m 3 = ^ _ p ' m 2 = m sC , wij = m 3 £ . (16) 

Thus the neutrino mass spectrum is fully determined. 

Table 1 shows the global-fit values of (^12^23^13) an d (5m 2 , Am 2 ) obtained by assuming 
CP conservation with cos<5 = ±1 [16]. We allow these parameters to independently vary in 
their 3a ranges: 

0.259 < sin 2 12 < 0.359 or 30.6° < 9 12 < 36.8° , 
0.340 < sin 2 23 < 0.640 or 35.7° < 9 23 < 53.1° , 
0.001 < sin 2 13 < 0.044 or 1.8° < 9 13 < 12.1° ; (17) 

and 

6.99 x 10" 5 eV 2 < 5m 2 < 8.18 x 10" 5 eV 2 , 

2.06 x 10" 3 eV 2 < |Am 2 | < 2.67 x 10~ 3 eV 2 . (18) 

In our numerical calculations we shall treat 5 as an unconstrained parameter. Note that the 
sign of Am 2 remains unknown: Am 2 > or Am 2 < corresponds to the normal or inverted 
mass hierarchy of three neutrinos. 



3 Two-zero Textures of M v 
3.1 Stability of texture zeros 

We first examine the stability of texture zeros of M v against the one-loop quantum corrections. 
To be explicit, we consider the unique dimension-5 Weinberg operator of massive Majorana 
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Table 1: The latest global-fit results of three neutrino mixing angles (#i 2) #23> #13) an d two 
neutrino mass-squared differences 5m 2 and Am 2 defined in Eq. (14). Here cos 5 — ±1, and 
the old reactor antineutrino fluxes have been assumed [16]. 



Parameter 


5m 2 (1(T 5 eV 2 ) 


Am 2 (10- 3 eV 2 ) 


#12 


#23 


#13 


Best fit 


7.58 


2.35 


33.6° 


40.4° 


8.3° 


la - range 


[7.32, 7.80] 


[2.26,2.47] 


[32.6°, 34.7°] 


[38.6°, 45.0°] 


[6.5°, 9.6°] 


2(7 range 


[7.16,7.99] 


[2.17,2.57] 


[31.6°, 35.8°] 


[36.9°, 50.8°] 


[5.1°, 10.9°] 


3<7 range 


[6.99,8.18] 


[2.06,2.67] 


[30.6°, 36.8°] 


[35.7°, 53.1°] 


[1.8°, 12.1°] 



neutrinos in an effective field theory after the heavy degrees of freedom are integrated out [18]: 



-d=5 

A 
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2«W 



-ah 



HH 



Toe 



+ h.c. , 



(19) 



where A is the cutoff scale, £ L denotes the left-handed lepton doublet, H = ia 2 H* with H 
being the standard- model Higgs doublet, and k stands for the effective neutrino coupling 
matrix. After spontaneous gauge symmetry breaking, H gains its vacuum expectation value 
(H) = v/y/2 with v ~ 246 GeV. We are then left with the effective Majorana mass matrix 
M v = kv 2 /2 for three light neutrinos from Eq. (19). If the dimension-5 Weinberg operator 
is obtained in the framework of the minimal supersymmetric standard model, one will be left 
with M v = k(v sin f3) 2 /2, where tan /3 denotes the ratio of the vacuum expectation values of two 
Higgs doublets. Eq. (19) or its supersymmetric counterpart can provide a way of generating 
tiny neutrino masses. There are a number of interesting possibilities of building renormalizable 
gauge models to realize the effective Weinberg mass operator, such as the well-known seesaw 
mechanisms at a superhigh energy scale A [19, 20, 21]. 

The running of M v from A to the electroweak scale /x ~ M z (or vice versa) is described 
by the renormalizat ion-group equations (RGEs) [22]. In the chosen flavor basis and at the 
one-loop level, M U (M Z ) and M J/ (A) are related to each other via 

(h 0\ (I e 0\ 



M V {M Z ) = I Q 





Vo 











MM) 





Vo 











(20) 



where the RGE evolution function 7 denotes the overall contribution from gauge and quark 
Yukawa couplings, and I a (for a — e, /i, r) stand for the contributions from charged-lepton 
Yukawa couplings [23]. Because of I e < I < I T as a consequence of m e <^ <^ m T , they 
can modify the texture of M v . In comparison, I ^ 1 only affects the overall mass scale of 
M v . Note, however, that the texture zeros of M v are stable against such quantum corrections 
induced by the one-loop RGEs. Taking Pattern A x of M u for example, we have 

/0 




V a 



(21) 



7 



at A, and thus 



M^(M Z ) 



( ° 




aIJ T \ 



K 



dI 2 T 



(22) 



at M z . This interesting feature implies that the important relations obtained in Eqs. (12) 
and (13) formally hold both at A and M z . In other words, if a seesaw or flavor symmetry 
model predicts a two-zero texture of M v at A, one may simply study its phenomenological 
consequences at M z by taking account of the same texture zeros. While the values of neutrino 
masses and flavor mixing parameters at M z turn out to be different from those at A, their 
correlations dictated by the texture zeros keep unchanged at any scale between A and M z . 



3.2 Permutation symmetry 

As pointed out in Refs. [10] — [14], the seven viable two-zero textures of M v can be classified 
into three distinct categories: (1) A 1 and A 2 ; (2) B l5 B 2 , B 3 and B 4 ; and (3) C. The 
phenomenological implications of those patterns in the same category have been found to be 
almost indistinguishable. Now we show that there exists a permutation symmetry between 
Patterns A x and A 2 , B x and B 2 , or B 3 and B 4 . This observation may help understand their 
similarities in model building and phenomenology. 

Let us take Patterns A 1 and A 2 in Eq. (2) for example. Note that the location of texture 
zeros in A 1 can be changed to that in A 2 by a permutation in the 2-3 rows and 2-3 columns. 
To be explicit, we define the elementary transformation matrix 



-^23 



/l 0\ 
1 

Vo l o/ 



(23) 



A A 

Then the Majorana neutrino mass matrix M v 2 can be constructed from M v 1 via 



M^=P 23 M^P? 3 . (24) 



If M v 1 is diagonalized by a unitary matrix V — UP like Eq. (8), then Eq. (24) tells us that 
M v 2 can be diagonalized by the unitary matrix V — UP with U = P 2 ?U . Parametrizing 
U in terms of three rotation angles (^12^23^13) an d one CP-violating phase 5, just like the 
parametrization of U in Eq. (9), we immediately obtain the relations 

#12 = ^12 > 013 = 013 > 023 = 2 ~~ ^23 > 5 = 5 -TV . (25) 

In addition, M^ 1 and Mjf" 2 have the same eigenvalues (for % = 1, 2, 3). The phenomenological 
predictions of Pattern A 2 can therefore be obtained from those of Pattern A-,^ by implementing 
the replacements in Eq. (25). It is easy to show that a similar permutation symmetry holds 
between Patterns B x and B 2 or between Patterns B 3 and B 4 . 
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3.3 Analytical approximations 

Thanks to the permutation symmetry discussed above, it is only necessary to study Patterns 
A l5 B-l, B 3 and C in detail. The analytical results for A 2 , B 2 and B 4 can be obtained, 
respectively, from those for A 1; B x and B 3 with the replacements 9 23 — > n/2 — 9 23 and 
5 — > 5 — 7r. In this subsection we explore the analytical relations among three neutrino masses 
(m 1 , m 2 , m 3 ), three flavor mixing angles {9 l2 , 9 23 , 9 l3 ) and three CP-violating phases (5, p, a) for 
each pattern of M u by means of Eqs. (11) — (16). The effective mass term of the neutrinoless 
double-beta decay, defined as (m) ee = \(M v ) ee \, will also be discussed. 

• Pattern A x with (Mj ee = (Mj e/1 = 0. With the help of Eq. (11), we obtain [11] 

^1 _ , s 13 | s 12 s 23 AS „ | 
A 3 c 13 \ c 12 c 23 / 

r = ~ s f (t 2 ? 6 '^ 5 ") ■ < 26) 

A 3 c 13 \*12 c 23 / 

Since sf 3 <C 1 is a good approximation, Eq. (26) immediately leads us to the results 



£ = — ~ tan 6 12 tan 9 23 sin 9 



m 3 
rrin 



13 5 



C = — w cot 6» 12 tan 6 23 sin 6» 13 ; (27) 

m 3 



and 



5 

ff W ■ (28) 

Without loss of generality, we choose three neutrino mixing angles to lie in the first 
quadrat and allow three CP-violating phases to vary in the ranges p G [— 7r/2, +7r/2], 
(j G [-7r/2,+7r/2] and 5 G [0,2tt]. Given 0.60 < tan# 12 < 0.75, 0.72 < tan# 23 < 1.3 and 
0.03 < sin# 13 < 0.21 at the 3<r level, Eq. (27) yields £ < ( < 1. So the normal neutrino 
mass hierarchy m 1 < m 2 < m 3 follows and Am 2 > holds. As indicated by Eq. (15), the 
Dirac CP-violating phase 5 can completely be fixed if three neutrino mixing angles and 
two neutrino mass-squared differences are known. In the lowest order of s 13 , however, 

= 5rrf_ _ 4 tan 2 9 23 sin 2 9 l3 
" Am 2 ~ sin2# 12 tan2# 12 ' 1 ' 

which is actually independent of 5. So we have to go beyond the leading order approxi- 
mation. To the next-to-leading order we obtain 



£ = — - ~ tan# 12 tan# 23 sin 9 13 \J 1 — 2 cot 9 12 cot 9 23 sin# 13 cos 5 , 

TTt I 

( = — - ~ cot # 12 tan 6* 23 sin # 13 Wl + 2 tan 6> 12 cot # 23 sin # 13 cos 5 ; (30) 

Wlr, 



9 



and 



R, 



C 2 



„ 4 tan 2 L sin 2 13 , . „ 

f = „:„o/i !—o/i (! + tan 2# 12 cot # 23 sin 13 cos 5) 



sin 29 12 tan 2$ 



12 



Finally we arrive at 
5 ~ 



cos 



tan 9. 



23 



/ sin2# 12 tan2# 19 .R 



12 J L i/ 



tan 26* 12 sin 13 y 4 tan 2 # 23 sin 9 13 



- 1 



(31) 



(32) 



Taking the best-fit values of the three neutrino mixing angles (i.e., 9 l2 = 33.6°, # 23 = 
40.4° and 9 l3 = 8.3°) together with those of two neutrino mass-squared differences (i.e., 
5m 2 = 7.58 x 10~ 5 eV 2 and Am 2 = 2.35 x 10~ 3 eV 2 ) from Table 1, one immediately 
obtains the neutrino mass spectrum 3 



TOo 



TO, 



TO 



4.8 x 10" 2 eV , 



m 3 cot 9 12 tan 9 23 sin 9 13 = 8.9 x 10 eV , 
m 3 tan 9 12 tan 9 23 sin 9 13 = 3.9 x 10~ 3 eV ; 



(33) 



and the CP-violating phases 5 ~ 64°, p ps 32° and <j —58°. Since (M^) ee = holds 
for Pattern A 1 of M u , the effective mass (m) ee of the neutrino less double-beta decay is 
definitely vanishing. 

• Pattern A 2 with {M u ) ee = (M u ) eT = 0. As pointed out in section 3.2, all the analytical 
results of Pattern A 2 can be obtained from those of Pattern A 1 with the replacements 
9 23 — > 7r/2 — # 23 and 5 — > 5 — ir. So it is straightforward to have 



f = — 

m 3 

m 2 
m 3 



tan 9 12 cot 9 23 sin 9 13 J 1 + 2 cot # 12 tan 6* 23 sin # 13 cos 5 , 



c 



cot # 12 cot 9 23 sin 9 13 \J\ — 2 tan 12 tan 9 23 sin 9 13 cos <5 ; 



and 



cos 



cot 0. 



23 



sin26> 12 tan 19 ^R 



12 ±L v 



tan 2# 12 sin 9 



'13 



4 cot 2 9 23 sin 2 9 13 



- 1 



(34) 



(35) 



In addition, 



5 7T 

9 ~ 2 2' 
47 ~ 2' 

and the neutrino mass spectrum turns out to be 



(36) 



TOo 



m 



TOi 



m 3 cot 12 cot # 23 sin 9 13 , 
m 3 tan 9 12 cot # 2 3 sm #13 • 



(37) 



3 Note again that such best-fit values have been obtained in the assumption of cos 6 = ±1 [16]. To give a 
rough estimate, however, we assume that they are essentially unchanged even for an arbitrary value of 8. 
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Similar to the case of Pattern A t , (m) ee = \(M u ) ee \ = holds in Pattern A 2 . Provided 
the same best-fit values of 5m 2 , Am 2 , 6 l2 , # 23 an d ^13 are taken, one can easily verify 
that cos 5 > 1 holds, implying a potential inconsistency of Pattern A 2 with current 
experimental data. When the uncertainties of those neutrino mixing parameters are 
taken into account (e.g., at the 3a level), however, we actually find no problem with 
Pattern A 2 (see the numerical analysis in section 4). 

Pattern B x with (M„) = (M u ) eT = 0. With the help of Eq. (11), we obtain 



e 2tS , 



■Sl2C 12 S 2 3(2c2 3 Sj i 3 S 23 C 13) c 23 S 13( S 12 s 23 e+t<5 + c 12 c 23 e ^) 
■^3 S 12 C 12 S 23 C 23 + ( S 12 — c 12) c 23 S 13 e * 5 + S 12 C 12 S 23 S 13 (1 + ^Lj)^ 

^2 _ S 12 C 12 S 23(^ C 23 g 13 ~ S 23 C 13) + C 23 g 13 ( C 12 g 23 e ~ H<5 + S 12 C 23 e ^) £ 2i5 ^gg^ 
^3 S 12 C 12 S 23 C 23 + ( S 12 — C 12) C 23 S 13 e * 5 + S 12 C 12 S 23 S 13 (1 + C 23) e2 * 5 

In the leading order approximation, 

i = — ~ tan2 ^23 , 



m 3 

C = ^ « tan 2 £ 23 ; (39) 
m 3 



and 



7T 

p Rj (J Rj S - - . (40) 



In the next-to-leading order approximation, we find 



m i m 2 _ _|_ 4 sin 13 cos 5 



m 3 m a sin 2# 12 sin 2# 23 



^ 2 sin 6 13 sin <5 . . 

^ °" sin 2# 12 tan 2# 23 tan 2 23 

The Dirac CP-violating phase 5 can be determined from 

^- 2S1 o! 13 |tan2fl 23 cos<5| . (42) 
sin zu^2 

Since 5m 2 > or equivalently m 2 > m x , we have cos 5 < from Eq. (41). Taking the 
best-fit values, one can obtain 5 RJ 91° and thus p rs a RJ 1°. Eq. (41) tells us that the 
difference between m 2 /m 3 and m 1 /m 3 is about 0.01, and that between p and <r is about 
4°. On the other hand, the neutrino mass spectrum is given by 

7.0 x 10~ 2 eV , 




m 3 tan 2 6 23 « 5.1 x 10~ 2 eV ; (43) 



ee ~ 



and the effective mass term of the neutrinoless double-beta decay turns out to be (m) 
m 3 tan 2 # 23 rs 5.1 x 10~ 2 eV. Note that we have input the best-fit value 8 23 = 40.4° [16], 
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so the normal neutrino mass hierarchy m 3 > m 2 > m l is allowed. If the 3a range of 9 23 is 
taken, however, both normal and inverted mass hierarchies are likely. When tan# 23 > 1, 
which is allowed at the 2a or 3a level [16], we have the inverted neutrino mass hierarchy 
with Am 2 < 0. In short, the analytical formulas for m 1 , m 2 and m 3 in Eq. (43) are 
valid no matter which mass hierarchy is taken, but the corresponding numerical results 
depend on the input value of 9 23 . 

Pattern B 2 with (M U ) TT = (M u ) = 0. By using the permutation symmetry discussed 
in section 3.2, we obtain 

f = — ~ COt2 9 23 , 



m 3 

C = ^ « cot 2 tf 23 ; (44) 



and 



7T 

p a « 5 — — (45) 

in the leading order approximation. Note that there is a period of n for p and a, so 
Eq. (45) is identical to Eq. (40) even if the replacement 5 — > 5 — n is made. In the 
next-to-leading order approximation, we have 

m 1 m 2 4 sin 6 l3 cos 5 

m 3 m 3 sin 29 12 sin 26> 23 

2 sin 9-, o sin 5 , ,„ s 

n-a^ ^ = . (46) 

sin 29 12 tan 29 23 cot 2 9 23 

The Dirac CP-violating phase is also determined by Eq. (42), because |tan2# 23 | keeps 
invariant under the replacement 9 23 — > n/2 — 9 23 . Now we have cos<5 > in view of the 
fact of m 2 > m 1 . So 5 m 89° and p « a m —1°. Eq. (46) tells us that the difference 
between m 2 /m 3 and m l /m 3 is about 0.01, and that between p and a is about 2°. The 
neutrino mass spectrum turns out to be 



m 



i ~ \l ~ 5.1 x 10~ 2 eV , 

3 Vl-cot 4 ^ 




7 23 

m 3 cot 2 23 « 7.0 x 10~ 2 eV . (47) 

In addition, (m) ee ~ m 3 cot 2 # 23 ~ 7.0 x 10~ 2 eV. Note that we have input the best-fit 
value 9 23 = 40.4° [16], so the inverted mass hierarchy m 2 > m 1 > m 3 appears. If the 
3a range of 9 23 is taken, however, both normal and inverted mass hierarchies are likely. 
Hence it is experimentally important to determine the deviation of 9 23 from 7r/4. 

Pattern B 3 with (M u )^ = {M v ) efX = 0. With the help of Eq. (11), we obtain 

^1 g 23 _ g 12 g 23 ~ C 12 C 23 g 13 e ^ ^iS 

^3 c 23 s 12 c 23 + c 12 s 23 s 13 e+ * 5 

^2 _ _£23 _ C 12 S 23 + S 12 C 23 S 13 e ^ c 2iS 

^3 c 23 c 12 c 23 — s 12 s 23' S 13 e+l<5 
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In the leading order approximation, 

m 
m 3 



771 

( = — « tan^ 23 ; (49) 



and 



7T 

p a 5 — — . (50) 



In the next-to-leading order approximation, we have 

m 1 m 2 4 tan 2 9 23 sin 9 13 cos 5 



m 3 m 3 sin 26* 12 sin 26> 23 

p-o-ftj + 13 . (51) 

^ sm20 12 tan 20 23 v y 

The Dirac CP-violating phase 5 can be determined from 

2 sin # 

R v « tan 2 # 23 Itan 2# 23 cos 51 . (52) 

" Sm20 12 23 1 23 I V ) 

The condition m 2 > m 1 leads to cos 5 > 0. Taking the best-fit values [16], we obtain 
5 ~ 89° and p ~ er —1° from Eq. (51). The difference between m 2 /m 3 and m 1 /m 3 
is about 0.008, and that between p and a is about 3°. On the other hand, the neutrino 
mass spectrum is given by 

7.0 x 10~ 2 eV , 
m 3 tan 2 9 23 « 5.1 x 10~ 2 eV . (53) 

The effective mass term of the neutrinoless double-beta decay turns out to be (m) ee « 
m 3 tan 2 # 23 5.1 x 10~ 2 eV. Note that the phenomenology of Pattern B 3 is essentially 
the same as that of Pattern B x except for the Dirac CP-violating phase 5. Note also that 
we have used the best-fit value 9 23 = 40.4° [16], so the normal neutrino mass hierarchy 
m 3 > m 2 > m 1 is allowed. If the 3a range of 9 23 is taken, however, both normal and 
inverted neutrino mass hierarchies are likely. 

Pattern B 4 with (M U ) TT = {M u ) eT = 0. By using the permutation symmetry discussed 
in section 3.2, we obtain 




23 ' 



and 



f = — « cot 2 9. 
m 3 

C = ^ « cot 2 £ 23 ; (54) 



5 - | (55) 
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in the leading order approximation. More accurately, we have 



m 



i. ± ~ _|. 



4 cot 2 # 2 3 sin 9 13 cos 5 



m 3 
p-o 



rrin 



sin 2# 19 sin 29 



'12 

2 sin 8 13 sin 5 



23 



sin 29 12 tan 29. 



23 



The Dirac CP-violating phase is given by 

2sin# 



R, 



— cot 2 ^3 | tan 2# 23 cos 8\ 



(56) 



(57) 



sin 2# 12 

We find cos 5 < from the fact of m 2 > m 1 . Hence 5 « 91° and p m a w 1°. Eq. (56) 
tells us that the difference between m 2 /m 3 and m l /m 3 is about 0.015, and that between 
p and a is about 3°. The neutrino mass spectrum turns out to be 




Am 2 



777.1 



COt 4 #23 
m 3 cot 2 #23 



5.1 x 10~ 2 eV , 
« 7.0 x 10~ 2 eV . 



(58) 



In addition, (m) ee ~ m 3 cot 2 # 2 3 ~ 7.0 x 10~ 2 eV. Note that we have used the best-fit 
value # 23 — 40.4° [16], so the inverted neutrino mass hierarchy m 2 > m 1 > m 3 appears. 
If the 3a range of # 2 3 is taken, however, both normal and inverted mass hierarchies are 
likely. It is obvious that a precise determination of 5 and the neutrino mass hierarchy is 
crucial to pin down one of the four patterns B 1; B 2 , B 3 and B 4 . All of them predict a 
nearly degenerate neutrino mass spectrum. 

Pattern C with (MJ = (Mj rr = 0. With the help of Eq. (11), we obtain 





C 12 C 13 




C 12( S 23 C 23) 


2 s 12 S 23 C 23 S 13 el<5 


A 3 


«13 


2Sl2 c 12 s 23 c 23 


~~ ( S 12 ~~ C 12)( S 23 


c 23J 6 13 e + Z6 12 c 12 6 23 c 23 6 13 e 


A 2 


s 12 c 13 




+ S 12( S 23 — C 23, 


— 2c l2 s 23 c 23 s 13 e l 


A 3 


s 13 


2 S 12 C 12 S 23 C 23 


~ ( S 12 — C 12)( S 23 


c 23J 6 13 e ^ ^*12 c 12*23 L 23*13 e 



AS 



(59) 



To the lowest order, 



£ = — 

m 3 



c = — 

m 3 



1 - 



2 cot 9 12 cos 5 



+ 



cot 2 9 



12 



\ tan 2#23 sin # 13 tan 2 26> 23 sin 2 # 13 



2 tan L cos 5 
1 + * . „ + 



tan 2 # 



12 



\ tan 2#23 sin 9 13 tan 2 2# 23 sin 2 # 13 



(60) 



Since m 2 > m 1 , Eq. (60) implies tan2# 23 cos<5 > and m 2 > m 3 . Furthermore, one can 
verify that tan 2# 12 tan 2# 23 sin # 13 cos S > 1 is required to guarantee m 2 > m 1 . Hence 
only the inverted mass hierarchy m 2 > m l > m 3 is allowed. The Dirac CP-violating 
phase 5 can be determined from 

2(1 + tan 9 12 tan #23) (t an 2# 12 tan 2# 2 3 sin # 13 cos 5 



1 



tan 9 l2 tan 9 23 tan 2# 12 (tan 9 12 + 2 tan 2# 2 3 sin # 13 cos S) 



(61) 



14 



Taking the best-fit values of three neutrino mixing angles and two neutrino mass-squared 
differences, we obtain 5 ~ 61° from Eq. (61). The Majorana CP-violating phases p and 
a turn out to be 



5 H — tan 
2 



5 tan 1 

2 



cot 9 12 sin 5 



tan 29 23 sin 9 l3 — cot 9 l2 cos S 



tan 9 l2 sin 5 



tan 29 23 sin 9 13 + tan 9 12 cos 5 



7T 

2 

7T 

2 



+13° 



-42° 



(62) 



Finally we obtain 



tan 2 29 23 cot 2 # 19 sin 2 9^ Am 2 



m 9 



12 1 



'13 z 



\ 1 + 2 cot # 12 tan 29 23 sin 9 13 cos 5 



4.3 x 10" 2 eV 



rrin 



2 tan # 12 cos 5 



\ tan 2^ 2 3 sin 9 13 



+ 



tern 2 9 



12 



tan 2 2#23 sin 2 9 13 



2 cot 9 12 cos 5 



+ 



cot 2 fl 



12 



^ tan 2^23 sin 9 13 tan 2 2# 2 3 sin 2 # 13 



6.6 x 10 -2 eV 



6.5 x 1(T 2 eV 



(63) 



together with 



1 



4 cot 2# 12 cos 5 



4 cot 2 29 



12 



tan 2 26U sin 2 9, 



4.1 x 10" 2 eV 



(64) 



t Jill '—9 23 0\ 3 i tui — l; 2.'i 1 1 ''13 

Note that both tan 2# 23 > with cos 5 > and tan 2# 23 < with cos 5 < are likely, if 
the 3a range of # 2 3 is taken into account. Nevertheless, the inverted mass hierarchy m 3 < 
m l < m 2 is expected in both cases. Going beyond the above analytical approximation, 
we find that the normal mass hierarchy m l < m 2 < m 3 is actually allowed if and only if 
#23 lies in the vicinity of 45° 4 . This point can be understood in a simple way. With the 
help of Eq. (59), one obtains 



Ai _ A2 

Aq Aq 



c 13 e 



2i8 



1 + s{ 3 e 



2i8 



-c 13 e 



z 2iS (l - sj 3 e 2lS ) 



(65) 



in the limit of 9 23 = 45°, implying the mass hierarchy m 1 = m 2 < m 3 for 9 13 7^ 0°. Hence 
a tiny deviation of 9 23 from 45° is required to lift the degeneracy of m 1 and m 2 and then 
reproduce the small but nonvanishing value of R u) such that m l < m 2 < m 3 comes out. 
Because the parameter space of this possibility is too small, we shall mainly concentrate 
on the inverted mass hierarchy in our subsequent numerical analysis. 

In summary, the analytical results for the three CP- violating phases (5, p, a), the two neutrino 
mass ratios m 1 /m 3 and m 2 /m 3 , the absolute neutrino mass m 3 and the effective neutrino mass 
(m) ee , predicted by seven two-zero textures of M u , are listed in Tables 2 — 4. If the best-fit 
values of neutrino mixing parameters [16] are taken, we find that only Pattern A 2 can be 
excluded. We expect that all the seven patterns of M v can survive current experimental tests 
at the 3cr level. A detailed numerical analysis will be done in section 4. 



J We thank the anonymous referee for pointing out this interesting possibility to us. 



15 



4 Numerical Analysis 



The analytical results obtained above show that the neutrino mass hierarchy is actually related 
to the flavor mixing angle 9 23 in Patterns B x , B 2 , B 3 and B 4 of M v . In particular, it depends 
on whether 9 23 > 45° or 9 23 < 45°. According to Table 1 [16], only 9 23 < 45° is allowed 
at the lcr level 5 . If a two-zero texture of M v can accommodate both normal and inverted 
mass hierarchies, we shall only concentrate on the one dictated by 9 23 < 45° in our numerical 
analysis, because the other possibility is just an opposite and trivial exercise. We have noticed 
that the Dirac CP-violating phase S should be close to 7r/2 or 37r/2 in Patterns T$ t 2 3 4 , and the 
differences between two Majorana CP-violating phases p and a in these two cases are distinct. 
For illustration, we shall only focus on the range of 5 around tt/2 in our numerical analysis of 
these four patterns. The explicit calculations are done in the following way: 

1. For each of the seven patterns of M v we generate a set of random numbers of (0 12 , 9 23) 9 13 ) 
and (5m 2 , Am 2 ) lying in their 3a ranges given by Eqs. (17) and (18) together with a 
random value of 5 in the range 5 G [0, 2ir]. 

2. Given the above random numbers, it is possible to calculate other physical parameters 
of M u , including three neutrino mass eigenvalues (m 1 ,m 2 ,m 3 ), two Majorana-type CP- 
violating phases (p,cr), the effective mass of the neutrinoless double-beta decay (m) ee 
and the Jarlskog invariant of leptonic CP violation J CP = Si2 c i2 s 23 c 23 s i3 c i3 s i n ^- To 
judge whether a pattern of M v is consistent with current experimental data or not, we 
require that the consistency conditions should be satisfied: (a) because of 5m 2 > 0, we 
require m 2 > m 1 or equivalently ( 2 — £ 2 > 0; (b) since only the neutrino mass hierarchies 
m 2 > m l > m 3 and m 3 > m 2 > m l are phenomenologically allowed, we further require 
(( 2 — 1)(£ 2 — 1) > (i.e., £ 2 > 1 and ( 2 > 1 correspond to the inverted mass hierarchy, 
whereas £ 2 < 1 and ( 2 < 1 stand for the normal mass hierarchy); (c) given the values of 
three neutrino mixing angles and two mass-squared differences, the Dirac CP-violating 
phase 5 is actually fixed by Eq. (15). Instead of solving 8, we require that Eq. (15) 
should be satisfied up to a reasonable degree of precision (e.g., 10~ 4 ). 

3. We consider all the points satisfying the consistency conditions, and then have a nine- 
dimensional parameter space spanned by nine quantities (9 12 , 9 23 , 9 13 , 5, p, a, m 1 ,m 2 , m 3 ). 
The low-energy observables such as J CP and (m) ee can accordingly be calculated. To 
present the final results in a simple and clear way, we restrict ourselves to the two- 
dimensional parameter space and set the x-axis to be the Dirac CP-violating phase 5. 
Therefore, what we actually show are the allowed ranges of relevant physical parameters 
of each pattern of M v . 

4. Corresponding to the allowed ranges of 9 12 , 9 23 and 9 13 changing with 5, their histograms 
are also plotted because they can signify the most probable values of three flavor mixing 

5 Note that 9 23 > 45° seems to be favored at the la level in the global analysis done by Schwetz et al [17]. 
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angles. The height of each histogram indicates the number of points in each bin. Since 
the future neutrino oscillation experiments will greatly improve the measurements of 
three neutrino mixing angles, such a presentation should be helpful in judging which 
two- zero pattern of M v is phenomenologically more favored. 

We stress that the strategy of our numerical analysis can also apply to other textures of the 
Majorana neutrino mass matrix M v . The two- zero patterns under discussion will serve as a 
good example to illustrate this strategy. 

Our numerical results are presented in Figs. 1 — 14. Comments and discussions follow. 

• Pattern A-^ Fig. 1 tells us that the neutrino mixing angles, in particular 9 l2 and 
#23, are actually insensitive to the Dirac CP- violating phase 5. This point can easily be 
understood with the help of Eq. (31) or Eq. (32), in which S is only loosely related to 
the ratio of two neutrino mass-squared differences R v due to the smallness of 9 13 . An 
interesting observation from the left panel of Fig. 1 is, that a relatively large value of 9 13 
(i.e., 9 13 ~ 7° • • • 8°), which is quite close to the best-fit value 6 l3 = 8.3° [16], is particularly 
favored. On the other hand, the neutrino mass spectrum is weakly hierarchical, as shown 
in Fig. 2. The dependence of m l and m 2 on 5 is ascribed to the next-to-leading order 
corrections given in Eq. (30). We have also illustrated the numerical prediction for J CP in 
Fig. 2. One can see the maximal value of J CP is at the percent level and should be able to 
lead to observable effects of CP violation in a variety of long-baseline neutrino oscillation 
experiments. Because 5 itself is essentially unconstrained by current experimental data 
at the 3a level, p and a turn out to be arbitrary as shown in Fig. 2, although their 
correlations with S are rather sharp. Finally we remark that Pattern A x of M v predicts 
(m) ee = for the neutrino less double-beta decay. 

• Pattern A 2 . As shown in Figs. 3 and 4, the phenomenological implications of Pattern 
A 2 of M v are essentially the same as those of Pattern A 1 . For instance, 9 13 w 6° ■ ■ -8° 
is favored and the effective mass term (m) ee is vanishing. Thus it is only necessary 
to emphasize their main difference. We have demonstrated a permutation symmetry 

A A 

between M v 1 and M v 2 in section 3.2 and found that the present best-fit values of 
neutrino mixing parameters (mainly 8 23 = 40.4° [16]) cannot coincide with Pattern A 2 . 
If the maximal mixing angle 9 23 = 45° were finally established, however, it would be 
almost impossible to distinguish between Patterns A 1 and A 2 in practice. 

• Pattern B x . Since m 2 > m 1 , we have cos 5 < (i.e., tc/2 < 5 < 3ir/2) in this case. 
In our numerical analysis we only focus on the range 5 G [tt/2,tt], because the range 
5 G [7T,37r/2] can similarly be analyzed. As shown in Figs. 5 and 6, only a very narrow 
region 5 G [0.507T, 0.567r] is phenomenologically allowed. This result obviously originates 
from Eq. (42), where | cos5| must be small enough to suppress the magnitude of R u . 
Furthermore, we only consider the normal neutrino mass hierarchy corresponding to 
9 23 < 45°. The latter seems to be favored by current data at the la level [16]. Fig. 5 
shows that 9 13 ~ 3° and 9 23 ~ 37° are more likely. A strong correlation between 9 23 and S 
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can also be understood from Eq. (42); namely, the maximal mixing 9 23 ps 45° requires the 
maximal CP- violating phase 5 ~ tt/2. In addition, a nearly degenerate mass spectrum 
as shown in Fig. 6 is predicted. There is a lower bound on (m) ee (i.e., (m) ee > 0.03 eV), 
which will be tested in the future experiments of the neutrinoless double-beta decay. 
Note that (m) ee « 0.2 eV can be achieved when 9 23 « 45° and 5 « 7r/2. The other three 
patterns of this category (i.e., B 2 , B 3 and B 4 ) have similar consequences, for which the 
numerical results have been given in Figs. 7 and 8, Figs. 9 and 10, and Figs. 11 and 12, 
respectively. But the details of these patterns, such as the allowed ranges of (p, a, 5) and 
(m 1 ,m 2 ,m 3 ), are somewhat different. 

• Pattern C. Fig. 13 shows no significant preference in the allowed ranges of three 
neutrino mixing angles. A strong correlation between 9 23 and S only appears when 5 is 
close to 7r/2. As shown in Fig. 14 and discussed in section 3.3, the inverted neutrino mass 
hierarchy is allowed in most parts of the parameter space (and a normal mass hierarchy 
is possible only when 9 23 is very close to 45° and 9 13 is nonvanishing) . Like Patterns 
Bi, 2,3,45 there is a lower bound on (m) ee in Pattern C (i.e., (m) ee > 0.02 eV), and its 
maximal value can saturate the present experimental upper bound (m) ee < 0.3 eV [24]. It 
should be noted that the sum of the three neutrino masses is subject to the cosmological 
bound J2 m i < 0-58 eV at the 95% confidence level [25], which has been derived from the 
seven-year WMAP data on the cosmic background radiation combined with the Baryon 
Acoustic Oscillations. Therefore, the possibility of 5 ~ tt/2 is essentially excluded at the 
same confidence level. Similar conclusions apply to Patterns B x 2 3 4 . 

Finally it is worth pointing out that we have also done a numerical analysis of the two-zero 
textures of M v by using the la and 2a values of 5m 2 , Am 2 , 9 12 , 9 23 and 9 13 (see Table 1 
and Ref. [16]). We find that all the seven patterns discussed above are compatible with 
current experimental data at the la or 2a level, although the corresponding parameter space 
is somewhat smaller. To be conservative, we take our numerical results obtained at the 3a 
level more seriously. 

5 Texture Zeros from Flavor Symmetries 

In general, the texture zeros of a Majorana neutrino mass matrix can be realized in various 
seesaw models with proper discrete flavor symmetries. It is possible to obtain the zeros in 
arbitrary entries of a fermion mass matrix by means of the Abelian symmetries (e.g., the 
cyclic group Z n [26]). To illustrate how to realize the two-zero textures of M v discussed above, 
we shall work in the type-II seesaw model, which extends the scalar sector of the standard 
model with one or more SU(2) L scalar triplets [20]. For N scalar triplets, the gauge-invariant 
Lagrangian relevant for neutrino masses reads 

- £ a = JEEW XZA^ L + h.c. , (66) 

J a,/3 
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where a and /3 run over e, \x and r, A ■ denotes the j-th triplet scalar field (for j — 1, 2, • • • , iV), 

and F A is the corresponding Yukawa coupling matrix. After the triplet scalar acquires its 

j 

vacuum expectation value (A ) = v A , the Majorana neutrino mass matrix is given by 

J 3 

M u = J2 Y a.v a ., (67) 

3 

where the smallness of v A is attributed to the largeness of the triplet scalar mass scale [20]. 

3 

In order to generate the texture zeros in M v and derive the seven viable patterns A 12 , 
B 1 2 3 4 and C, we follow the spirit of Ref. [27] and impose the Z n symmetry on the Lagrangian 
in Eq. (66). The unique generator of the cyclic group Z n is w = e* 27r// ™, which produces all the 
group elements Z n = {1, w, w 2 , ■ ■ ■ , w 11 " 1 }. Now it is straightforward to specify the number of 
scalar triplets N (i.e., j = 1, 2, • • • , N), the order of the cyclic group n and the representations 
of £ a and A ■ under the symmetry group. 

• N = 3 and n = 6 for A 12 and B 3 4 . In this case we have to introduce three scalar triplets 
and the Z 6 symmetry group, for which the generator is w = e m ^. The representations 
of lepton doublets £ ah (for a = e,fi,r) and the scalar triplets A^ (for j = 1,2,3) are 
assigned as follows. 

1. For Pattern A x : 

4l^ 2 4l, ^L^l^L, £ Th ^W% L , 

A 1 ->■ 1 A 1 , A 2 ->■ zu 3 A 2 , A 3 ->■ zuA 3 . (68) 

Given the above representations of lepton doublets, three scalar triplets are needed 
to enforce the nonzero elements in the neutrino mass matrix: A 1 for (M v ) and 
(M U ) TT , A 2 for (MJ MT and A 3 for (Mj er . 

2. For Pattern A 2 : 

£ eL ^ w 2 £ eL , e^-nfa, i TL ^w\ L , 

A 1 ->■ 1 A ± , A 2 ->■ zu 3 A 2 , A 3 ->■ ro 4 A 3 . (69) 

Note that Eq. (69) differs from Eq. (68) only in the assignment for the triplet A 3 . 
Such a difference originates from the fact that (M v ) efI = and {M v ) eT ^ hold for 
Pattern A l5 while {M u ) eT = and (M u ) efM ^ hold for Pattern A 2 . It is worthwhile 
to point out that the assignments in Eq. (69) are by no means unique. As we have 
discussed in section 3.2, there exists a permutation symmetry between A x and A 2 . 
Therefore, we may exchange the representations of £^ lL and £ rL in Eq. (68) but 
preserve those of scalar triplets to obtain A 2 from A 1 . 

3. For Pattern B 3 : 

£ eL ^l£ eL , £, L ^ru\ L , £ rL ^w% L , 

A 1 ->■ 1 A ± , A 2 ->■ w 3 A 2 , A 3 ->■ zuA 3 . (70) 
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4. For Pattern B 4 : 

e eL ^w% L , ^41^, e rL ^w\ L , 

A 1 -> 1 A 1 , A 2 ro 3 A 2 , A 3 ->• ro 4 A 3 . (71) 

In addition to the permutation symmetry P 23 discussed in section 3.2, we notice 
that the location of texture zeros in Pattern B 3 can be obtained from that in 
Pattern A x by a permutation in the 1-2 rows and 1-2 columns. Similarly there 
exists a 1-3 permutation symmetry between Pattern B 4 and Pattern A 2 . Although 
these symmetries cannot lead to simple relations between any two of the neutrino 
mixing parameters, they are instructive for the assignments of lepton doublets. For 
instance, Eq. (70) and Eq. (71) can be derived from Eq. (68) and Eq. (69) by 
exchanging the representations of £ eL and £^ L and the representations of £ eL and £ rL , 
respectively. The assignments of scalar triplets should not be changed. 

• N = 2 and n = 3 for B 12 . In this case we introduce only two scalar triplets and the 
Z 3 symmetry group, for which the generator is w = e 2l7r//3 . The assignments of lepton 
doublets £ aL (for a — e, /i, r) and the scalar triplets A ■ (for j = 1, 2) are as follows. 

1. For Pattern B x : 

'eL-U'eL, ~+ , U ™Xl , \ ~+ 1 \ , A, ^ ^A 2 . (72) 

2. For Pattern B 2 : 

£ eL ^l£ eL , * t l-^tl, A^lAn A 2 ^u7 2 A 2 . (73) 

Note that the assignment in Eq. (73) is slightly different from that in Ref. [27], 
where A 2 — > wA 2 is taken and the representations of lepton doublets are the same 
as in Eq. (72). Here we have implemented the 2-3 permutation symmetry between 
Pattern B x and Pattern B 2 . 

• N = 3 and n = 4 for Pattern C. In this case we have to introduce three scalar triplets 
and the Z 4 symmetry group, for which the generator is w = e l7r//2 . The assignments of 
lepton doublets £ aL (for a = e, //, r) and the scalar triplets A ■ (for j = 1, 2, 3) are 

Aj-^lAn A 2 ->■ zu 3 A 2 , A 3 ->■ zuA 3 . (74) 

Thus all the seven two-zero patterns of M v can be obtained in this simple symmetry scheme. 

In all cases we have taken the right-handed charged-lepton singlets E aR to transform in the 
same way as the left-handed lepton doublets £ aL and taken the standard-model Higgs doublet 
H to be in the trivial representation. Hence the charged-lepton mass matrix M l is diagonal, 
as we have chosen from the beginning. The two- zero textures of M v can also be realized in the 
seesaw models with three right-handed neutrinos, several Higgs singlets, doublets and triplets, 
by imposing either Abelian or non-Abelian discrete flavor symmetries [26, 27, 28]. 
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6 Summary 



In view of the latest T2K and MINOS neutrino oscillation data which hint at a relatively large 
value of 9 13 , we have performed a systematic study of the Majorana neutrino mass matrix M v 
with two independent texture zeros. It turns out that seven patterns (i.e., A 12 , B 1j2i3j4 and 
C) can survive current experimental tests at the 3a level, although they are also compatible 
with the data at the la or 2a level. The following is a brief summary of our main observations: 

• Given the values of three flavor mixing angles (9 12 , 9 23 , 9 l3 ) and two neutrino mass-squared 
differences (5m 2 , Am 2 ), it is in principle possible to fully determine three CP- violating 
phases (8, p, a) and three neutrino masses (m 1 ,m 2 ,m 3 ). The analytical formulas for the 
latter have been derived and listed in Tables 2 — 4. 

• By making the analytical approximations and taking the best-fit values of neutrino mixing 
parameters [16], we find that only Pattern A 2 can be excluded. We have numerically 
confirmed that all the seven patterns of M v (i.e., A 12 , B 1 2 34 and C) are compatible 
with current neutrino oscillation data at the la level, but our numerical results have been 
presented only at the more conservative 3a level. 

• Figs. 1 — 14 show the main numerical results of our systematic analysis. Some interesting 
points should be emphasized. (1) Both A 1 and A 2 favor a relatively large 9 13 (e.g., 
9 l3 ~ 8°), B x 234 prefer a relatively small 9 13 (e.g., 9 l3 ~ 3°), and C shows no significant 
preference for the magnitude of 9 13 . (2) The Dirac CP- violating phase 5 obtained from 
B 1 2 34 lies in a narrow range around tt/2 or 3ir/2, and 5 = 7f/2 is strongly correlated 
with 9 23 — 7r/4. (3) For 5 — > n/2 and 9 23 — > 7r/4, the predictions for neutrino masses m i 
and the effective neutrino mass {m) ee may run into contradiction with their upper bounds 
set by the cosmological observations and the neutrinoless double-beta decay experiments. 
(4) The size of J CP may reach the percent level and thus appreciable leptonic CP violation 
is possible to show up in the future long-baseline neutrino oscillation experiments. 

In addition we have shown that the texture zeros of the Majorana neutrino mass matrix M v 
are stable against the one- loop quantum corrections 6 , and pointed out that there exists a 
permutation symmetry between A x and A 2 , B : and B 2 or B 3 and B 4 . In the type-II seesaw 
model with two or three scalar triplets we have illustrated how to realize two- zero textures of 
M v by using the Z n flavor symmetry. 

The ongoing and upcoming neutrino oscillation experiments are expected to measure the 
neutrino mixing parameters, in particular the smallest mixing angle 9 l3 , the deviation of 9 23 
from 7r/4 and the Dirac CP-violating phase 5. The sensitivity of future cosmological observa- 
tions to the sum of neutrino masses m i and the sensitivity of the neutrinoless double-beta 
decay experiments to the effective mass term (m) ee will probably reach ~ 0.05 eV in the near 

6 In contrast, the texture zeros of the Dirac neutrino mass matrix are essentially sensitive to quantum 
corrections like those of quark mass matrices [29]. 
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future. We therefore expect that some patterns of the two-zero textures of M v might be ex- 
cluded or only marginally allowed by tomorrow's data, and those capable of surviving should 
shed light on the underlying flavor structures of massive neutrinos. 
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Table 2: Seven viable patterns of the neutrino mass matrix M v with two texture zeros, and 
their predictions for three CP- violating phases (5, p, a). 



Pattern Texture of M„ 



CP-violating phases 



/ x \ 

x x 

\ X X X / 



COS 



tan 9. 



23 



sin 26 12 tan 29 l2 R u 



tan 29 12 sin 9 A s \ 4 tan 2 6* 9S sin 2 6* 



- 1 



'13 



23 1 



'13 



OH 

H 2 ' 2 2 



/ x \ 

XXX 

V x x/ 



cos 



cot 9 



23 



tan2# 12 sin# 



5 
2 



'13 
7T 

2 



sin 2# 12 tan 29 12 R h 
4 cot 2 # 23 sin 2 9 l3 

5_ 

2 



a 



fx x \ 

x x 
V x x/ 



cos 



sin 2^, 



2 sin # 13 | tan 2^23 1 



5 ~2 ' 



2 sin 13 sin 5 



sin 2^ 1 2 tan 20 23 tan 2 fl 23 



fx x \ 

x x 
\x x 0/ 



cos 



sin2# 19 i? 



12 1L u 



2 sin# 13 | tan 29 



A 71 

S ~2> 



23 I 



2 sin # 13 sin 5 



sin 26 1 2 tan 2fl 23 cot 9 23 



/ x x \ 
x 

\ X X X / 



COS 



A 71 

5 ~2 ' 



sin 2# 12 cot 2 # 23 K 



2sin# 13 | tan 2^23 

p — (J Rj 4 



2 sin # 13 sin 5 



sin 2flj 2 tan 29 23 



fx x \ 

xxx 
V x / 



cos 



A 71 
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sin 2# 12 tan 2 # 23 K 
2 sin # 13 | tan 2^3 1 

p — (J Rj -t 



2 sin # 13 sin 5 
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C 



2(1 + tan^tantf 
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tan 2 # 12 tan 2# 12 tan 9 23 R U 
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fx X X \ 

x x 
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[1 + (1 — i?^) tan 6* 12 tan 6* 23 ] tan 2# 12 tan 2# 23 sin 6* 13 



r 1 1 

5 + - tan 1 



5 tan 1 

2 



cot # 12 sin S 



tan 2#23 sin 9-, , — cot 9 cos 5 



'13 

tan 9 12 sin <5 



12 1 



tan 26> 23 sin 6> 13 + tan 9 12 cos 5 



7T 

2 

7T 
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Table 3: Seven viable patterns of the neutrino mass matrix M v with two texture zeros, and 
their predictions for two neutrino mass ratios £ = m 1 /m 3 and ( = m 2 /m 3 . 



Pattern Texture of M„ 



Neutrino mass ratios 




£ pa tan 9 12 tan 9 23 sin 9 13 , C ~ c °t #12 tan # 2 3 sin # 13 



/ x 
x x 
V x 



0\ 
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12 
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Table 4: Seven viable patterns of the neutrino mass matrix M v with two texture zeros, and their 
predictions for the absolute neutrino mass m 3 and the effective mass terms of the neutrinoless 
double-beta decay (m) ee . 



Pattern Texture of M v The scales of neutrino masses 





/ x \ 

x x 

\ X X X / 
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Figure 1: Pattern A 1 of M u : allowed ranges of flavor mixing angles (6' 12 , 6 , 23 , 6' 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 2: Pattern A x of M v : allowed ranges of the neutrino masses (m 1 ,m 2 ,m 3 ), the Jarlskog 
invariant J CP and the Majorana CP-violating phases (p, a) versus the Dirac CP-violating phase 
S at the 3<T level. 
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Figure 3: Pattern A 2 of M u : allowed ranges of flavor mixing angles (6' 12 , 6 , 23 , 6' 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 4: Pattern A 2 of M u : allowed ranges of the neutrino masses (m 1 ,m 2 ,m 3 ), the Jarlskog 
invariant J CP and the Majorana CP-violating phases (p, a) versus the Dirac CP-violating phase 
S at the 3<T level. 
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Figure 5: Pattern B x of M v : allowed ranges of flavor mixing angles (6> 12 , 9 23 , 9 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 6: Pattern B x of M v : allowed ranges of the neutrino masses (m 1 ,m 2 ,m 3 ), the Jarlskog 
invariant J CP and the Majorana CP-violating phases (p, a) versus the Dirac CP-violating phase 
5 at the 3a level. 
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Figure 7: Pattern B 2 of M^: allowed ranges of flavor mixing angles (6> 12 , 9 23 , 9 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 8: Pattern B 2 of M v : allowed ranges of the neutrino masses (m 1 ,m 2 ,m 3 ), the Jarlskog 
invariant J CP and the Majorana CP-violating phases (p, a) versus the Dirac CP-violating phase 
5 at the 3a level. 
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Figure 9: Pattern B 3 of M v : allowed ranges of flavor mixing angles (6> 12 , 9 23 , 9 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 10: Pattern B 3 of M v : allowed ranges of the neutrino masses (m 1 , m 2 , m 3 ), the Jarlskog 
invariant J CP and the Majorana CP- violating phases (p, <r) versus the Dirac CP- violating phase 
S at the 3(7 level. 
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Figure 11: Pattern B 4 of M v : allowed ranges of flavor mixing angles (#i 2) #23> $13) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 12: Pattern B 4 of M v : allowed ranges of the neutrino masses (m 1 , m 2 , m 3 ), the Jarlskog 
invariant J CP and the Majorana CP- violating phases (p, a) versus the Dirac CP- violating phase 
S at the 3a level. 
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Figure 13: Pattern C of M u : allowed ranges of flavor mixing angles (6> 12 , 9 23 , 9 13 ) versus the 
Dirac CP-violating phase 5 at the 3a level, where the probability distribution of three angles 
are shown in the left panel. 
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Figure 14: Pattern C of M v : allowed ranges of the neutrino masses (m 1 ,m 2 ,m 3 ), the Jarlskog 
invariant J CP and the Majorana CP-violating phases (p, a) versus the Dirac CP-violating phase 
S at the 3a level. 
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